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INTRODUCTION 


Physical  phenomena  and  mechanisms  central  to  the  degradation  of  epoxy 
resins  and  epoxy  resin  composites  have  been  extensively  examined  over  the 
two  years  of  this  Grant*  The  experimental  and  theoretical  methods  used 
have  led  to  the  Identification  of  factors  that  need  to  be  considered  when 
substituting  composites  for  metals  as  primary  structural  members*  The 
programme  began  with  attention  directed  towards  the  prevention  and 
prediction  of  loss  of  mechanical  strength,  specific  attention  being  paid 
to  the  further  development  and  application  of  optical  methods  for 
monitoring  loss  of  ability  to  transfer  shear  stress  at  Interfaces  between 
fibre  and  matrix  materials.  As  the  study  has  progressed.  It  has  become 
Increasingly  apparent  that  direct  optical  methods  offer  reliable  means  for 
assessment  of  the  Importance  of  changes  occurring  at,  or  very  close  to, 
the  interface. 

One  technique  which  shows  considerable  promise  Is  that  in  which-  the  load 
transferred  from  resin  to  fibre  Is  measured  photoelastlcally.  A  microscope 
with  rectified  optics  has  Improved  our  facility  for  making  these 
measurements  and,  by  paying  careful  attention  to  specimen  preparation, 
recent  results  have  permitted  comparisons  to  be  made  between  fibres  of 
different  composition  and  different  surface  condition.  The  method  used 
here  Is  the  only  direct  way  that  such  comparisons  can  be  made,  and  further 
development  of  the  technique  to  enable  the  experiment  to  be  controlled  by 
means  of  a  microcomputer  could  pay  enormous  dividends. 

Ultrasonic  photoelasticity  investigations  have  continued  In  an  attempt  to 
acquire  more  precise  data  from  property  changes  attributable  to  water 
uptake.  By  the  very  nature  of  Its  modes  of  propagation,  ultrasound 
samples  a  very  much  larger  region  of  material  than  that  which  can  be 
sensibly  identified  as  "the  Interface"  and,  for  this  reason,  the  technique 
Is  more  sensitive  to  macroscopic  changes  such  as  matrix  plasticlsatlon 
than  It  Is  to  localised  changes  such  as  the  destruction  of  bonds  at 
fibre/matrix  Interfaces.  Efforts  to  raise  the  frequency  and  sharpen  the 
ultrasonic  fringe  system  are  unlikely  to  contribute  to  the  Interfacial 
dialogue. 

Work  in  this  area  is  very  relevant  to  non-destructive  evaluation  (NDE)  of 
composites  In  the  field.  In  respect  of  this,  analysis  of  data  collected 
from  several  transducers  giving  piezoelectric  output,  present  challenging 
problems  In  both  acquisition  and  analysis.  Stress  fields  arising  from 
differential  contraction  between  fibre  and  resin  materials  and 
modifications  of  these  stress  fields  that  result  from  resin  swelling 
during  water  uptake  are  now  analysed  by  finite  element  method.  The  finite 
element  method  (fern)  due  to  Adams  and  co-workers  (1970)  treats  the  purely 
elastic  fibre  and  elastoplastlc  matrix  In  generalized  plane  strain  and 
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provides  for  a  pseudo-three-dimensional  stress  analysis*  The  method 
adopted  here  is  that  developed  by  Alexander  and  Turner  (1975)  for 
axlsymmetrlcal  flow  and  has  been  modified  in  order  to  accommodate  plane 
strain*  The  interface  is  represented  by  a  third  phase*  The  analysis  so 
far  attempted  parallels  the  methods  reported  by  Adams* 

Two  new  optical  methods  have  been  explored  in  order  to  secure  experimental 
measurement  of  internal  stresses  associated  with  Inhomogeneous  resin 
shrinkage  and  swelling.  The  first  is  an  oblique  incidence  polarising 
microscope  method  which  has  been  used  to  uniquely  determine  the  principal 
stresses  transverse  to  the  fibre  direction  in  an  as-cured  uniaxial 
composite*  The  second  method  exploits  the  optical  Interference  pattern 
created  between  a  resin  layer  (actually,  between  a  thin  microscope  cover 
slip  in  contact  with  the  resin  layer)  and  an  optical  flat  in  order  to 
measure  the  resin  swelling.  Using  thin  plate  elasticity  theory,  the  resin 
displacements  so  measured  have  been  converted  to  stresses* 

Large  scale  programs  have  been  run  on  the  GEC  40/85  minicomputer  and 
smaller  programs,  one  of  which  receives  data  typed  directly  into  an  on¬ 
line  computer  console,  use  the  Cromemco  System  3  micro.  The  latter 
computer,  now  situated  in  our  composites  laboratory,  is  proving 
Invaluable.  Plans  for  the  interfacing  of  the  micro  with  the  piezoelectric 
transducer  for  sampling  a  wave-form  very  rapidly  (~16  MHz)  are  at  an 
advanced  stage. 


2*  DEBONDING  OF  FIBRES  IN  AN  EPOXY  RESIN  MATRIX 


2.1  POINCARE 

The  stresses  generated  in  a  glass  fibre-epoxy  resin  composite  during  resin 
cure  and  during  subsequent  cooling  to  room  temperature  are  small  in 
magnitude  and  complex  in  distribution.  Their  analysis  using  polarized 
light  is  simplified  by  representing  the  state  of  polarization,  and  the 
change  in  this  state  as  the  light  passes  through  the  optical  system,  by 
means  of  the  PoincarS  sphere  construction,  Poincar4(1892).  A  single  point 
on  the  surface  of  a  sphere  represents  the  polarization  of  electromagnetic 
oscillation  at  any  point  in  the  system.  The  changes  in  this  polarization 
are  represented  by  the  motion  of  the  representative  point  about  the 
surface  of  the  sphere  in  accordance  with  simple  rules.  In  a  polarized 
light  wave  the  end  point  of  the  electric  vector  JJ  describes  an  elliptical 
helix.  This  helix  projects  on  to  a  plane  parallel  to  the  wave-front  as  an 
ellipse  tAiose  eccentricity  is  described  by  the  angle  $  ■  tan-1  (b/a), 
where  a  and  b  are  the  major  and  minor  axes.  The  orientation  of  the  ellipse 
with  respect  to  some  arbitrary  Cartesian  basis  (i,  j)  normal  to  the  optic 
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axis  is  specified  by  the  angle  a  (Figure  1).  The  point  P  representing  the 
state  of  polarization  described  by  the  angles  a  and  0  is  plotted  on  the 
surface  of  a  unit  sphere  at  a  latitude  of  2$  relative  to  the  equator  and  a 
longitude  of  2a  relative  to  some  arbitrary  meridian  (Figure  2). 
Conventionally,  points  in  the  upper  hemisphere  represent  electromagnetic 
radiation  In  which  the  £  vector  rotates  antl-clockwise,  and  those  In  the 
lower,  clockwise* 

This  mode  of  representation  has  several  consequences  that  are  Immediately 
obvious*  All  points  lying  on  the  equator  represent  linearly  polarized 
light;  the  two  poles  (North  and  South)  represent  circularly  polarized 
light  rotating  anti-clockwise  and  clockwise  respectively;  points  on  the 
same  parallel  (0  ■  constant)  represent  ellipses  of  the  same  elliptlclty 
and  different  orientations  and  conversely  points  on  the  same  meridian  (a  ■ 
constant)  represent  ellipses  of  varying  form  but  the  same  orientation; 
two  points  at  the  opposite  ends  of  the  same  diameter  represent  elliptical 
helices  of  opposite  handedness  -  both  ellipses  have  the  same  form  but  the 
major  axis  of  one  ellipse  coincides  with  the  minor  of  the  other* 

This  last  proposition  has  an  Important  special  case  triien  the  points  lie  on 
the  equator,  and  therefore  represent  two  waves  linearly  polarized  at  right 
angles*  These  waves  can  generate  any  elliptically  polarized  wave  and  the 
Polncar£  sphere  can  be  used  to  determine  the  required  phase  and 
amplitude  relationships  between  them*  In  Figure  2,  I  and  J  represent  the 
two  linearly  polarized  waves  and  P  represents  the  elliptically  polarized 
wave  to  be  generated.  Conversely,  P  represents  a  polarized  wave, 
elliptlclty  angle  0,  trtiich  is  to  be  resolved  into  two  orthogonal  linearly 
polarized  components  in  the  directions  .1  and  ,j,  a  clockwise  and  90-a 
anti-clockwise  from  the  major  axis  (Figure  1).  In  both  cases  the 
construction  Is  the  same* 

The  15  vector  of  an  elliptically  polarised  e-m  wave  as  a  function  of  time 
Is 


E^(t)  "  a  cos  cot  %  +  b  sin  cot  £  (1) 

where  £  and  ?  are  unit  vectors  parallel  to  the  major  and  minor  axes  whose 
lengths  are  a  and  b  respectively  (See  Figure  1).  Using  axes  JL  and  j 
Inclined  at  a  to  8  and  E  can  be  expressed 

15 (t)  *  A  cos  (u)t-9j)_l  +  B  cos  (cot-02 )«L  (2) 

where  6j  and  @2  are  the  phase  angles  and  A  and  B  the  amplitudes  of  the  two 
waves*  To  compare  these  two  expressions  we  express 
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(3) 


55  •  cos  a  j.  +  sin  a  j 
£  ■  -sin  a  _i  +  cos  a  j 
Substituting  (3)  into  (1) 


X  "  [a  cos  wt  cos  a  -  b  sin  wt  sin  a]^ 

4-  [a  cos  wt  sin  o  +  b  sin  wt  cos  a]J^ 
Expanding  (2)  and  comparing  terms  with  equation  (4)  yields 
a  cos  a  ■  A  cos  0^  b  sin  a  -  -A  sin  0^ 


(4) 


(5) 


a  sin  a  -  B  cos  0 


2  b  cos  a  ■  B  sin  ©2 


(6) 


By  eliminating  0^  and  ©2  the  following  equations  are  easily  obtained 
A2  4  B2  «  a2  +  b2 
A2  -  B2  -  (a2  -  b2)  cos  2a 
Putting  B/A  -  tan  <J>  these  two  equations  can  be  reduced  to 
cos  2<|)  ■  cos  2(3  cos  2a 
In  the  spherical  triangle  IPQ 

cob  IOP  ■  cos  2(3  cos  2a 

hence  24>  ■  IOP  which  determines  the  relative  amplitudes  of  the  two  waves. 
From  equations  (5)  and  (6)  can  be  obtained 

2ab  -  2AB  sin(0o  -  0,)  (a2-b2)  tan  a  «  2AB  008(00-0.) 

TT77  T775  2  1  7TT? 


a  +  b*1  A  4-  B 


A  -  B 


which  on  rearranging,  become 

sin  2(3  -  sin  24>  8in  (@2  -  0j)  tan  2a  -  tan  2<\>  cos  (02  -  0j) 
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In  the  spherical  triangle  IPQ 

sin  20  »  sin  2<l>  sin  0  tan  2a  •  tan  2tp  cos  0 

where  0  is  the  angle  between  the  equatorial  plane  and  the  plane  containing 
I,  J  and  P. 

Hence  0  ■  ©2  -  0j,  indicating  that  the  wave  represented  by  J  lags  by  0 
behind  the  wave  I* 

On  entering  a  birefringent  crystal,  polarized  light  can  be  considered  as 
split  into  two  components,  parallel  to  the  fast  and  slow  axes  of  the 
crystal  in  the  plane  normal  to  the  optic  axis.  Their  phase  and  amplitude 
relationships  are  given  by  the  angles  0  and  <\>  described  in  the  previous 
paragraph.  As  the  two  waves  pass  through  the  crystal  at  different  speeds 
their  relative  amplitudes  remain  the  same  but  the  phase  difference  between 
them  changes  by  an  angle  (60)  proportional  to  the  thickness  (e)  of  crystal 
traversed  and  also  to  the  difference  in  the  refractive  index  of  the 
crystal  for  light  waves  polarized  in  the  fast  and  slow  directions  (An). 
In  the  Poincarfe  sphere  representation,  4>  remains  constant  and  9  changes  by 
60  -  2fcAne/X  radians.  Assuming  I  and  J  to  be  the  fast  and  slow  directions 
respectively  (Figure  3),  the  position  of  P,  which  represents  the  incident 
light,  is  rotated  by  60  anti-clockwise  about  the  axis  IJ  as  indicated. 

The  Brace-Kohler  compensator  is  used  to  measure  small  optical 
retardations,  particularly  in  connection  with  photo-elastic  experiments. 
It  consists  of  a  slice  of  birefringent  crystal  which  introduces  a  known 
retardation  RQ  between  light  waves  polarized  in  the  fast  and  slow 
directions  of  the  crystal.  It  is  usually  calibrated  for  a  given 
wavelength  of  light  and  the  retardation  is  normally  quoted  in  nanometres 
of  path  difference.  A  typical  experimental  arrangement  of  crossed  polars, 
compensator  and  specimen  is  shown  in  Figure  4.  Initially  the  compensator 
is  aligned  with  its  fast  and  slow  axes  parallel  to  the  axes  of  the  crossed 
polars  and  the  fast  axis  of  the  specimen  is  set  at  45°  to  the  polars.  The 
compensator  is  then  rotated  through  an  angle  <J>  until  the  intensity  of 
light  passing  through  the  whole  system  is  a  minimum;  the  retardation  R  due 
to  the  specimen  is  then  given  by  the  following  expression: 

R  "  RQ  sin  24< 

The  Poincar4  sphere  analysis  of  the  above  optical  system  is  shown  in 
Figure  5.  P  and  A  represent  the  directions  of  the  polarizer  and  analyser 
respectively.  As  discussed  in  the  previous  section,  the  point  P  denoting 
light  emerging  from  the  polarizer  is  rotated  by  the  compensater  by  an 
angle  0  equal  to  2nRQ/X  about  an  axis  CC'  lying  in  the  equatorial  plane  at 
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an  angle  24>  to  the  line  PA;  hence  P  moves  to  P'.  The  specimen  rotates  the 
point  P'  through  an  angle  p  equal  to  2xR/\  about  the  axis  SS' 
(perpendicular  to  AP)  to  a  point  P"  such  that  POP"  is  a  minimum,  the 
condition  for  minimum  light  intensity  transmitted  through  the  analyser. 
Hence  P"  must  lie  on  the  equatorial  circle.  If  we  assume  that  OP  -  1,  then 
it  is  easily  shown  that 

cot  9  sin  P  •  cos  0  sin  24> 

from  which  we  obtain 

sin  P  «*  sin  0  sin  2<J> 

For  small  retardations,  the  following  approximation 
R  *  Rc  sin  2<l> 
holds,  as  before. 


2.2  DISCUSSION  OF  THE  EXPERIMENTS 

When  an  epoxy  resin  cures,  and  as  it  cools  from  the  cure  temperature,  it 
shrinks  and  generates  stresses  in  and  around  inclusions  such  as  totally 
embedded  short  fibres.  The  stress  fields  are  cylindrically  symmetric  with 
radial  compressive  stresses  in  the  resin  transmitted  normally  across  the 
fibre/resin  interface.  The  resin  adjacent  to  the  fibre  also  experiences 
hoop  tensile  stresses  and  a  tensile  stress  parallel  to  the  fibre  which  is 
under  axial  and  radial  compression.  Figure  6  shows  the  pattern  of  stress- 
induced  birefringence  in  and  around  an  'S'  glass  fibre  drawn  in  ammonia 
and  totally  embedded  in  a  slab  of  Ciba  Geigy  MY  750  epoxy  resin;  the  fibre 
axis  is  aligned  parallel  to  one  of  the  crossed  polars.  The  stress  within 
the  fibre  includes  a  component  which  arises  from  shear  stress  transferred 
at  the  interface  and  which  builds  up  with  distance  from  each  fibre  end. 
The  existence  of  axial  shear  stresses  in  the  resin  adjacent  to  the  fibre 
ends  is  demonstrated  by  the  lobes  of  light  extending  a  third  of  the  length 
of  the  fibre  from  both  ends.  The  series  of  photographs  in  figure  7  taken 
between  crossed  polars  at  different  specimen  orientations  (in  steps  of  5°) 
gives  an  indication  of  how  the  purely  tensile  stresses  in  the  resin 
parallel  to  the  fibre  at  its  centre  gradually  acquire  a  shear  component 
towards  the  fibre  ends.  The  load  transfer  index  is  a  measure  of  the  build 
up  of  compressive  stresses  in  the  fibre  due  to  interfacial  shear  stresses. 
It  is  expected  that  as  the  interfacial  bonds  break  down  the  compressive 
stress  build  up  will  be  reduced  until,  when  total  debonding  has  occurred, 
there  will  be  no  variation  in  fibre  compression  along  its  length.  Figure  8 
is  a  schematic  representation  of  the  stresses  in  the  fibre  and  resin. 
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fhe  load  tfartsfer  Index  (LTI)  for  a  given  fibre  18  defined  as  the 
difference  between  the  corrected  retardations  measured  in  the  middle  of 
the  fibre  and  the  average  of  the  corrected  retardations  measured  one  fibre 
diameter  in  from  both  ends.  The  corrected  retardation  is  the  difference 
between  the  retardation  measured  through  the  fibre  and  the  average  of  the 
retardations  measured  adjacent  to,  and  one  fibre  diameter  from,  the  fibre. 
All  compensation  measurements  are  taken  with  the  fibre  at  45°  to  the 
crossed  polar  axes  in  accordance  with  the  procedure  for  measuring 
retardations  using  the  Brace-Ktthler  compensator  described  in  the 
discussion  of  the  Poincar&  sphere.  During  the  last  nine  months  a  re¬ 
appraisal  of  the  measurement  of  the  LTI  from  a  theoretical  and 
experimental  standpoint  has  led  to  a  fuller  appreciation  of  the  wealth  of 
Information  that  can  be  obtained  from  retardation  measurements  through  and 
adjacent  to  embedded  fibres.  Since  a  number  of  unquantlf lable  factors 
affect  the  absolute  values  of  retardation  and  load  transfer  index,  the 
major  Importance  of  the  LTI  lies  in  its  variation  and  gradual  decline 
after  successive  periods  of  Immersion  in  warm  water.  However  the 
Polncar&  sphere  construction  has  increased  our  understanding  of  the 
significance  of  the  individual  measurements. 

The  specimens  studied  were  all  cast  in  silicon  rubber  moulds.  A  suitable 
quantity  of  glass  fibre  was  cut  into  1-2  mm  lengths,  washed  in  acetone, 
dried  in  an  oven  and  added  to  the  liquid  epoxy  resin  prior  to  vigorous 
mixing  and  subsequent  de-gassing  under  vacuum.  The  moulds  were  carefully 
cleaned  and  dried,  and  warmed  to  approximately  50°C  to  aid  in  filling  the 
mould  to  the  desired  level  with  the  viscous  resin.  After  gelling  for  2 
hours  at  100°C  the  oven  temperature  was  raised  slowly  to  150°C  and  the 
specimens  were  cured  for  a  further  4  hours  before  slow  cooling  in  the 
oven.  Despite  these  precautions,  all  the  specimens  possessed  intrinsic 
birefringence  with  the  fast  axis  approximately  parallel  to  the  long  axis 
contraction  between  resin  and  mould.  The  magnitude  of  this  background 
birefringence  was  at  least  double  the  greatest  stress-induced 
birefringence  arising  from  the  resin  shrinking  onto  a  fibre,  hence  the 
resultant  pattern  of  isoclinlcs  on  and  around  a  fibre  only  demonstrated 
the  expected  cylindrical  symmetry  when  the  fibre  was  approximately  aligned 
with  the  fast  or  slow  axes  of  the  background  birefringence.  When  this 
condition  obtained  the  fast  axes  in  the  resins  either  side  of  the  fibre 
end  were  misoriented  from  the  fibre  axis  by  approximately  equal  but 
opposite  angles. 

Using  the  Poincarh  sphere  construction  it  can  easily  be  shown  that  if  the 
slow  and  fast  axes  of  the  area  of  specimen  being  studied  are  not  at  45°  to 
the  crossed  polars,  the  angle  of  compensation  is  reduced.  However,  if  the 
misorientation  angles,  although  of  opposite  sign,  are  equal  in  magnitude 
(as  described  above),  the  corrections  to  the  compensation  readings  are  the 


same*  This  was  confirmed  by  the  following  measurements  taken  from  a 
commercial  'E'  glass  fibre  aligned  at  6°  to  the  specimen  axis,  at  which 
orientation  the  retardations  wither  side  of  a  fibre  end  were  19*5  nm  and 
17*0  nm  respectively.  Rotating  the  specimen  19.3°  clockwise  brought  the 
first  area  Into  correct  alignment  and  a  reading  of  25  nm  was  recorded.  A 
rotation  of  21.2°  anti-clockwise  changed  the  second  retardation  reading  to 
23.25  nm,  hence  the  two  readings  were  changed  approximately  in  proportion. 
It  is  interesting  to  note  that  if  the  Inverse  cosines  of  the  two  ratios 
between  uncorrected  and  corrected  retardations  are  taken,  the  results  are 
38.7°  and  43°  respectively  which  are  twice  the  misalignment  angles,  as 
predicted  by  the  Poincar4  sphere  construction. 

Photographed  between  crossed-polars,  the  fibre  in  Figure  9  was  surrounded 
by  a  symmetrical  pattern  of  isoclinics  similar  to  Figure  6  but  the  fibre 
Itself  was  not  completely  dark.  Rotating  the  specimen  4°  clockwise 
rendered  the  fibre  dark  but  destroyed  the  cylindrical  symmetry  of  the 
resin  pattern.  This  is  an  example  of  slight  misalignment  between  the 
stress-induced  and  background  birefringence.  In  more  extreme  cases,  it 
was  not  possible  to  orient  the  fibre  so  that  the  isocllnics  displayed 
cylindrical  symmetry.  The  isoclinics  were  always  zero  order  but  because 
of  the  background  birefringence  they  did  not  necessarily  indicate  regions 
of  the  specimen  in  which  the  principal  stress  directions  were  aligned  with 
the  crossed  polars. 

The  fibres  studied  were  of  four  types  -  commercial  ' E '  and  'S'  glass,  'S' 
glass  drawn  in  a  dry  ammonia  atmosphere  and  'E'  glass  drawn  in  a 
tetraf luoromethane  (CF^)  atmosphere  (these  two  laboratory  glasses  were 
supplied  by  David  Martin);  in  the  main  experiment  two  of  each  type  of 
fibre  were  studied.  The  commercial  fibres  were  all  10  microns  in  diameter 
but  for  the  laboratory  fibres  the  'E'  glass  fibres  were  6  microns  and  the 
'S'  were  13  microns.  Each  measurement  of  LTI  entailed  a  minimum  of  nine 
compensation  readings,  three  at  each  end  and  three  in  the  middle  of  the 
fibre.  A  corrected  retardation  was  also  obtained  approximately  midway 
between  the  middle  and  one  end.  In  addition  the  background  compensation 
in  the  resin  away  from  the  Influence  of  any  fibre  induced  stresses  was 
recorded  together  with  the  '*balance"  compensation  at  which  the  intensity 
of  light  passing  through  the  fibre  centre  was  the  same  as  that  passing 
through  adjacent  resin. 

All  the  retardations  are  theoretically  calculable  as  integrals  along  the 
light  path  and  correspond  to  more  or  less  complex  curves  on  the  Poincarfe 
sphere.  The  retardation  of  light  passing  through  the  fibre  has  components 
due  to  the  differences  between  (1)  the  hoop  and  axial  stresses  (both 
tensile)  in  the  resin,  and  (il)  the  radial  end  axial  compressive  stresses 
in  the  fibre.  Retardation  in  the  resin  ne/t  to  the  fibre  centre  is  a 
relatively  simple  case  as  the  axial  tension  and  the  resultant  of  hoop  and 
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radial  stresses  are  always  aligned  with  the  crossed  polars*  At  the  fibre 
ends*  the  situation  Is  complicated  by  the  shear  stresses  in  the  resin 
which  result  In  the  principal  stress  directions  being  misaligned  with  the 
fibre,  as  discussed  above.  The  magnitudes  of  these  various  stresses  are 
determined  Initially  by  the  fibre  geometry,  the  difference  In  thermal 
expansion  coefficients  between  resin  and  glass  and  the  position  of  the 
fibre  within  the  slab  of  epoxy  resin  (the  fibres  tend  to  settle  between  5 
and  20  microns  from  the  bottom  of  the  specimen).  By  symmetry  the 
interfacial  shear  stress  at  the  centre  of  the  fibre  Is  zero.  Mismatch 
between  fibre  and  resin  as  a  result  of  resin  shrinkage  Increases  with 
distance  from  the  fibre  centre,  x,  producing  shear  stresses  which, 
according  to  Cox  (1952),  are  proportional  to  slnh(x).  The  radial 
compressive  and  hoop  tensile  stresses  are,  according  to  simple  theories, 
equal  In  magnitude  and  both  fall  off  as  1/r  where  r  Is  the  radial  distance 
from  the  fibre  axis.  Their  magnitudes  were  expected  to  depend  on  fibre 
thickness,  being  greatest  for  the  13  micron  laboratory  'S'  glass  fibres 
and  least  for  the  6  micron  laboratory  'E'  glass.  This  was  borne  out  by 
experimental  observation  in  that  the  retardation  In  the  resin  next  to  the 
fibre  centre  differed  from  the  background  retardation  by  11.5  nm  for  a 
laboratory  'S'  glass  fibre,  5.85  nm  for  commercial  'S'  glass  and  1.5  nm 
for  laboratory  'E'  glass.  The  variation  of  retardation  with  radial 
distance  from  the  axes  of  the  laboratory  'S'  glass  fibre  is  shown  In 
Figure  10  and  shows  the  expected  1/r  dependence. 

The  proximity  of  the  fibres  to  the  specimen  surface  allows  stress 
relaxation  at  the  surface  to  occur  reducing  the  magnitude  of  the  stresses 
generated.  It  Is  noticeable,  when  considering  the  raw  data  that  the 
thinner  the  fibre,  the  greater  the  variation  In  retardation  along  the 
length  of  the  fibre  and  the  less  the  variation  of  the  retardation  In  the 
resin.  This  indicates  a  greater  stress  variation  In  the  thinner  fibres 
due  to  its  smaller  cross-sectional  area.  In  addition  It  Is  known  that 
exposing  the  epoxy  resin  to  water  alters  its  stress  optic  coefficient  and 
causes  swelling.  Work  is  in  progress  to  Investigate  these  phenomena.  At 
this  stage  some  Idea  of  the  magnitude  of  these  effects  can  be  gleaned  from 
the  variation  In  background  birefringence  which  changed  by  a  maximum  of 
5.85  nm  for  any  one  specimen  during  the  complete  run  and  the  typical 
change  was  of  the  order  of  1  nm  between  Immersions.  1  nm  of  retardation 
In  polarised  light  passing  through  50  microns  of  epoxy  resin  indicates  a 
principal  stress  difference  of  3.5  x  10^N/m2  assuming  a  stress  optic 
coefficient  of  56  Brewsters. 
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2.3 


SUMMARY  OF  RESULTS 


A  decision  was  taken  early  In  the  second  year  of  work  to  carry  out  the 
water  Immersion  tests  at  60°C.  It  was  found  that  If  the  specimens  were 
boiled,  the  debonding  process  occurred  in  a  catastrophic  fashion  with 
little  significant  variation  of  load  transfer  index  before  the  sudden  and 
total  loss  of  interface  integrity.  As  an  example  of  this  behaviour,  the 
SYMVU  plot  (Figure  11)  of  corrected  retardation  at  various  positions  along 
an  'E'  glass  fibre  drawn  in  CF^  against  immersion  time  in  boiling  water 
shows  that  after  15  minutes  the  fibre  'lost  contact'  with  the  resin  and 
there  was  as  a  result  no  variation  of  stress  along  the  length  of  the 
fibre.  A  similar  test  on  commercial  'E'  glass  indicated  total  loss  of 
load  transfer  after  13  minutes,  an  insignificant  variation  considering  the 
individual  immersions  were  of  longer  duration.  It  was  expected  that  at 
60°C  the  processes  of  deterioration  would  progress  fast  enough  to  be 
observable  over  a  realistic  time  span.  In  the  light  of  data  collected 
since  that  decision,  it  is  appropriate  for  further  work  to  be  carried  out 
at  75°C.  Above  80°C  the  diffusion  coefficient  of  water  in  MY  750  departs 
from  the  normal  Boltzmann  variation  with  temperature.  At  60°C  an 
overnight  soaking  frequently  produced  little  or  no  effect  on  the  fibres. 

Figure  12  summarises  the  data  obtained  from  the  specimens  during  170  hours 
of  immersion  in  water  at  60°C.  The  large  variations  in  load  transfer 
index  which  occurred  in  all  cases  during  the  first  1500  minutes  is  not 
wholly  attributable  to  experimental  error  which  was  +2  nm  in  the  LTI. 
The  general  downward  trend  of  each  graph  is  overlaid  with  fluctuations 
which  arise  from  the  change  in  stress  optic  coefficient  and  resin  swelling 
as  described  above.  The  thinner  glass  fibres  ('E'  in  CF^)  appeared  to 
have  less  load  transferred  to  them  and  the  'S'  glass  in  ammonia  fibres 
initially  had  higher  values  of  LTI  than  the  commercial  fibres.  However 
during  immersion,  several  bubbles  developed  on  the  surface  of  the  two  'S'- 
in-ammonla  fibres,  indicating  total  debonding  particularly  at  the  fibre 
ends.  Figure  13  is  an  optical  micrograph  of  one  of  these  fibres  after 
6400  minutes  immersion.  The  load  transfer  index  reflected  this 
deterioration,  falling  to  3.4  nm  for  the  fibre  in  figure  13  and  1.1  nm  for 
the  other.  The  other  fibres  all  experienced  the  initial  slight  drop  in 
LTI  but  during  the  final  80%  of  the  immersion  time  there  was  no 
significant  change,  because  time  was  running  short  the  specimens  were 
placed  in  boiling  water  for  1  minute  and  then  left  in  the  water  as  it 
slowly  cooled.  The  effect  on  this  treatment  was  in  some  cases 
considerable  as  can  be  seen  from  the  micrographs  in  figures  14,  15  and  16. 
In  particular  the  'E'  glass  fibre  drawn  in  CF^  (figure  14)  debonded  almost 
completely  towards  the  ends.  The  commercial  'S'  glass  fibres  were  least 
affected  by  the  boiling  treatment.  The  bubble  visible  at  the  end  of  the 
fibre  in  Figure  15  existed  prior  to  boiling  which  served  to  increase  its 
size.  It  seems  likely  that  the  boiling  caused  small  invisible  bubbles  of 
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water  to  form  in  regions  of  high  water  content  and  coalesce  into  the  large 
visible  bubbles  seen  in  the  micrographs*  The  test  pointed  up  the 
unsuitability  of  boiling  as  the  normal  test  mode.  The  commercial  ' E ' 
glass  fibre  in  figure  16  was  more  affected  than  commercial  'S'  glass,  with 
the  level  of  debonding  rising  to  about  50%  at  the  fibre  end.  It  was 
noticeable  that  the  bubbles  formed  preferentially  in  the  regions  of 
highest  shear  stress  at  the  ends  of  the  fibres.  This  may  have  been  a 
factor  in  the  rapid  debonding  of  the  'S'  glass  drawn  in  ammonia  since 
experiments  showed  that  these  fibres  had  the  highest  levels  of  interfacial 
shear  stress. 

The  results  of  the  initial  tests  carried  out  during  the  first  year  of  work 
are  summarised  in  Table  1.  These  were  in  the  nature  of  exploratory  tests 
and  indicate  the  need  for  the  more  sophisticated  approach  adopted  in  the 
second  year. 


2.4  DISCUSSION  OF  THE  RESULTS 

The  research  has  shown  that  the  load  transfer  experiment  is  a  useful  and 
necessary  simplification  of  a  complex  physical  process.  The  subtraction 
of  the  retardation  in  resin  adjacent  to  the  fibre  does  not  afford  a 
perfect  correction  to  the  retardation  through  both  resin  and  fibre  but 
results  have  shown  that  the  load  transfer  index  measured  in  this  way 
provides  a  simple  way  of  assessing  the  load  bearing  capabilities  of  the 
fibre  resin  interface.  Points  germane  to  the  experimental  technique  are 
discussed  below,  and  this  is  followed  by  a  critical  evaluation  of  the 
results,  emphasising  their  relevance  to  assessing  the  suitability  of 
various  fibre  chemistries  for  reinforcing  epoxy  resins. 

It  was  found  vital  that  all  fibres  studied  underwent  the  same  regime  of 
preparation  and  immersion.  In  practice  it  was  a  simple  matter  to  place 
all  the  specimens  in  the  same  water  bath  thus  fulfilling  the  second 
criteria.  However,  the  preparation  of  the  specimens  posed  more  logistical 
problems  since,  from  a  batch  of  six  epoxy  resin  specimens  containing  short 
fibres,  there  were  on  average  only  one  or  two  fibres  per  specimen  worthy 
of  study.  In  section  2.3  some  of  the  factors  that  might  influence  the 
value  of  the  LTI,  when  no  changes  in  the  properties  of  the  Interface  had 
occurred,  were  discussed.  These  included  variation  in  the  resin  stress 
optical  coefficient  with  water  content,  the  stresses  produced  by 
inhomogeneous  resin  swelling  due  to  water  Ingress  and  the  position  of  the 
fibre  relative  to  the  surface  of  the  specimen.  Rotating  the  mould  in  the 
oven  during  gel  would  create  a  more  random  distribution  of  the  fibres  in 
the  specimen  but  would  also  lead  to  problems  due  to  resin  flow  before  and 
during  gel,  and  also  due  to  fibres  lying  at  an  angle  to  the  specimen 
surface.  This  would  complicate  the  interpretation  of  the  measured 
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RESIN 

FIBRE 

WATER 

TIME  TAKEN 

TEMPERATURE 

TO  DEBOND 

MY  750 

Carbon  'A' 

22°C 

no  debonding  in  81  days 

Carbon  'A' 

100°C 

“  9  hours 

Carbon  'H' 

22°C 

no  debonding  in  81  days 

Carbon  'H' 

100°C 

15  minutes 

'E'  glass  (NH3  drawn) 

22°C 

40-53  days 

'E'  glass  (NH3  drawn) 

22°C 

40-53  days 

coupled* 

'E'  glass  (NH^  drawn) 

100°C 

<  1  1/2  hours 

'E'  glass  (NH3  drawn) 

100°C 

1  3/4  hours 

coupled* 

'S'  glass  (NH3  drawn) 

22°C 

36-49  days 

'S'  glass  (NH3  drawn) 

22°C 

36-49  days 

coupled* 

'S'  glass  (NH3  drawn) 

100°C 

22  hours 

'S'  glass  (NH3  drawn) 

100°C 

22  hours 

coupled* 

SP  250 

S2  glass 

22°C 

25  days 

S2  glass 

80°C 

1-5  days 

MY  720 

'S'  glass  (NH3  drawn) 

22°C 

9  days 

'E'  glass  (NH3  drawn) 

22°C 

9  days 

456  (S2)  glass 

22°C 

8  days 

P263A  (S2)  glass 

22°C 

8  days 

801  (E)  glass 

22°C 

not  debonded  in  8  days 

H  +  S  904  (SI)  glass 

22°C 

not  debonded  in  8  days 

S  glass  (NH3  drawn) 

80°C 

<  20  1/2  hours 

456  (S2)  glass 

80°C 

1-2  hours 

P263A  (S2)  glass 

80°C 

1  hour 

801  (E)  glass 

80°C 

1-2  hours 

H  +  S  904  (SI)  glass 

80°C 

1-4  days 

*  Fibres  coated  with  silane  coupling  agent,  A187,  manufactured  by  Union 
Carbide 


TABLE  1*  Loss  of  load  transfer  by  various  composite  specimens 

during  exposure  to  water 


12 


retardations*  The  ideal  specimen  would  contain  fibres  of  the  same  uniform 
shape  and  size  embedded  at  equal  distances  from  the  surface  in  otherwise 
stress-free  resin*  However,  knowledge  of  the  changes  induced  by 
departures  from  this  ideal,  together  with  Informed  observation  of  the 
visible  effects  of  water  attack,  allow  realistic  comparisons  to  be  made* 

The  'S'  glass  fibres  drawn  in  the  laboratory  contain,  or  are  coated  with, 
water  soluble  constituents,  dissolution  of  which  material  creates  pressure 
pockets  at  the  fibre/matrix  interface  after  170  hours  of  exposure  of  the 
composites  to  60®C  water*  The  fact  that  premature  failure  is  not  observed 
with  laboratory  drawn  'E'  glass  fibres  (l*e*  failure  which  is  premature 
compared  with  the  behaviour  of  commercially  drawn  'E'  glass  fibres)  is 
attributed  to  the  differences  in  chemical  nature  of  the  gaseous  atmosphere 
into  which  drawing  takes  place*  It  is  understood  that  Professor  Martin 
used  a  tetraf luoromethane  (CF^)  atmosphere  for  his  'E'  glass  fibres  and  an 
ammonia  (NH^)  atmosphere  for  his  'S'  glass  fibres*  Ammonium  salts  are 
readily  soluble  in  water,  saturated  aqueous  solutions  of  which  generate 
osmotic  pressures  having  magnitudes  in  excess  of  100  bars  at  the  boiling 
point  and  probably  also  at  60°C.  The  adhesive  strength  of  glass /wet  epoxy 
resin  joints  at  60®C  is  expected  to  be  of  the  same  order  of  magnitude* 
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3 


ULTRASONIC  PHOTOELASTICITY 


The  photoelastic  method  for  Imaging  ultrasound  was  developed  for 
investigating  the  material  property  changes  at  or  near  glass  interfaces. 
Turner  and  Ashbee  (1978).  Use  of  a  light  emitting  diode  controlled  by 
TTL,  Andrews  and  Wallis  (1977),  gave  pulse  widths  of  50  ns,  and  by 
measuring  velocities  of  the  longitudinal  and  shear  waves,  it  is  possible 
to  observe  the  Influences  of  interfacial  healing,  continuing  cure  of 
heated  specimens,  plasticisatlon  on  water  uptake  and  debonding  on 
weathering,  Farrar,  Turner,  and  Ashbee  (1979).  Such  data  give  information 
about  the  bulk  material  and  details  of  the  properties  of  the  Interfacial 
bond  are  not  available.  Experimental  difficulties  encountered  include 

(a)  rapid  attenuation  of  the  ultrasound  signal  and 

(b)  poor  resolution  arising  from  the  frequency  which,  at  2  MHz,  is 
somewhat  low  for  collecting  data.  It  was  therefore  decided  that  the 
interface  might  be  more  easily  studied  in  specimens  of  single  crystal 
quartz  and  epoxy  resin.  Changes  in  mechanical  damping  characteristics 
during  water  uptake  should  be  revealed  by  using  such  resin  coated  quartz 
piezoelectric  cyrstals  to  drive  the  transducer  in  visualization 
experiments.  The  large  amount  of  data  contained  in  photographs  is  not 
easily  extracted  from  the  wave  trace  and  the  analysis  must  await  the 
computer  interface  which  will  enable  the  changes  in  damping 
characteristics  to  be  analysed  completely. 

3.1  QUARTZ  OSCILLATOR 

The  state  of  localised  stress  near  the  interfacial  region  that  is 
introduced  during  curing  precludes  the  direct  observation  of  ultrasonic 
fringes  in  the  resin.  Successful  application  of  the  technique  is 
restricted  to  observation  of  fringes,  in  an  adjacent  stress  free 
visualizing  block,  i.e.  remote  from  the  specimen.  In  order  to  investigate 
Si02/epoxy  resin  interfaces  more  directly,  a  series  of  experiments  have 
been  carried  out  where  the  PZT5A  piezoelectric  discs  are  replaced  by 
crystalline  quartz.  X-cut  slabs,  2  mm  thick,  were  sawn  from  a  synthetic 
quartz  single  crystal  having  a  cross-sectional  area  of  about  2  cm.  The 
slabs  were  then  polished,  cleaned  and  coated  with  1  mm  of  Ciba-Geigy  MY750 
epoxy  resin.  An  electric  field  up  to  800v.  was  applied  across  the 
epoxy /quartz /epoxy  with  aluminium  foil  electrodes,  and  then  shorted  using 
a  thyrister.  The  resonant  mechanical  vibration,  observed  on  an 
oscilloscope,  begins  with  an  amplitude  of  15v.  and  decays  in  about  15 
cycles.  The  measured  frequency  is  typically  39  MHz  which  is  more  than  an 
order  of  magnitude  higher  than  the  expected  frequency  of  1.44  MHz.  The 
complexity  of  the  signal  is  made  up  of  a  series  of  modes  of  vibration  of 
the  quartz  together  with  the  reflections  from  the  epoxy  coating  with  the 
major  component  '1.5  MHz,  expected  from  the  thyrister  voltage  cut-off 
rate,  'lost'  in  the  total  waveform.  The  interfacial  bonding  on  boiling  the 
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specimens  in  distilled  water  modifies  the  frequency  trace  by  changing  the 
relative  pulse  heights.  The  energy  emanating  from  the  quatz  slabs  so  far 
investigated  is  not  sufficient  to  drive  a  transducer  for  ultrasound 
visualisation,  but  the  modified  trace  is  there  for  analysis.  It  was 
decided  to  analyse  the  signal  by  means  of  the  computer.  Sampling  the  form 
of  the  wave  by  means  of  a  fast  A/D  converter  is  regarded  as  the  most 
effective  method  and  the  technique  to  do  this  has  been  designed. 

Work  in  this  area  was  suspended  pending  the  filing  of  a  patent 
application,  however,  owing  to  unintentional  public  release  of  our  final 
technical  report  under  Grant  DA-ERO-76-G-068  (R-D  No  2160),  the  patent 
application  has  not  been  filed. 


3.2  COMPUTER  MODEL 

The  work  on  ultrasonic  photoelasticity  has  contributed  less  to  Interfacial 
physics  for  one  principal  phenomenon:  stress  waves  can  never  be  observed 
using  the  visualisation  techniques.  Turner  &  Ashbee  (1978),  in  a  stress 
gradient.  It  was  nevertheless  quickly  realised  that  the  ultrasound  could 
pick  up  useful  data,  in  a  mode  converting  reflection,  and  pass  on  into  a 
stress  free  zone  where  the  fringes  could  be  observed.  Experimentation 
with  composites  acoustically  coupled  to  a  visualising  quartz  block  was 
found  to  be  effective  although  attenuation  particularly  of  the  transverse 
wave  is  significant.  Steel  specimens  similarly  coupled  exhibit  less 
attenuation  since  the  electron  and  crystal  dislocation  scatter  in  metals 
is  smaller  than  that  produced  in  two  phased  composites.  Work  on  the  steel 
specimens  proceeded  along  NDE  lines  and  has  been  reported  elsewhere, 
Kitson,  Low  and  Turner  (1979).  The  computer  model  originally  designed  for 
composite  specimens  and  described  by  Turner  &  Ashbee  (1978),  has  since 
been  applied  to  the  non-destructive  application  with  developments  to 
improve  the  graphics  on  our  own  computer  equipment  with  the  interface 
between  the  specimen  and  visualising  block  taken  into  account. 
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4 


STRESS  FIELDS 


4.1 


CURING  STRESSES 


Differential  thermal  contraction  between  fibre  and  matrix  materials  during 
cooling  from  the  cure  temperature  is  not  accommodated  by  resin  flow  and 
the  state  of  self-stress  attributable  to  this  source  has  been  calculated 
for  simple  fibre  geometries  by  Adams  and  co-workers  (1979)  &  (1977).  They 
assume  that  the  contraction  is  uniformly  distributed  and  in  glass 
fibre/epoxy  resin  systems,  their  finite  element  model  predicts  that  the 
radial  principal  stress  is  tensile  along  the  line  of  centres  between 
adjacent  fibres  and  compressive  elsewhere,  and  that  its  magnitude  at  the 
fibre/resin  interface  is  in  the  range  100  MPa  to  -80  MPa. 

The  aim  of  this  aspect  of  the  present  research  is  to  measure  two  of  the 
principal  components  of  the  self-stress  field  of  a  glass  fibre/epoxy  resin 
composite  in  the  as-cured  state. 

In  order  to  reduce  the  problem  to  one  of  two-dimensional  elasticity,  thin 
transverse  sections  were  cut  from  a  long  unidirectional  composite. 
Accommodation  of  shrinkage  by  way  of  rigid  body  displacements  in  the 
transverse  plane  were  prevented  by  casting  the  composite  inside  a  thick 
walled  glass  tube.  The  unrelaxed  radial  and  hoop  principal  stresses 
within  any  slice  subsequently  removed  from  the  specimen  were  then  amenable 
to  measurement  by  the  oblique  incidence  photoelastic  method  described  by 
Durelli  and  Riley  (1965). 

The  two-dimensional  stress-optic  law  for  a  beam  of  light  that  enters  the 
surface  of  a  thin  transverse  section  at  normal  incidence  is 


R 


-  o2) 


(7) 


where  Oj  and  o2  are  the  principal  stresses  in  the  transverse  plane 
\  »  wavelength  of  the  light 
C  ■  stress  optical  coefficient 
t  ■  thickness  of  transverse  section 
In  terms  of  the  fringe  order  n,  equation  (1)  becomes 


n 


Tmax 


(8) 


tjnax  *  maximum  shear  stress  in  the  transverse  plane 
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Equation  (8)  can  also  be  written  as 

*-*' 

where  F  -  X 

2Ct 

is  the  model  stress  fringe  value. 

The  model  stress  fringe  value  is  the  change  in  maximum  shear  stress 
necessary  to  produce  a  change  of  unity  in  the  fringe  order.  The  quantity 


(9) 

(10) 


Ft  -  \m  f  (11) 

Is  independent  of  the  thickness  of  the  model  and  is  called  the  material 
stress  fringe  value.  f  is  the  change  in  maximum  shear  stress  required  to 
produce  a  change  of  unity  in  fringe  order  in  a  model  of  unit  thickness. 
Substituting  equation  (5)  into  equation  (2), 


n 


(a,  -  a2)  t 
27 


Let  the  specimen  be  rotated  about  principal  axis  1,  through  an  angle  0. 
The  path  length  (tQ)  of  the  light  becomes 


t_  -  t 
°  cos© 


(12) 


and  the  principal  stresses  transform  to 


Let 


Oj  ;  c?2  cosZ6 

al  *  ill’  a2  ” 

2f  It 


The  corresponding  fringe  orders  identified  by  nQ^  are 

n0<”  *  (Qt  -  o9  cosZ9)t  "  (aT  -  a?  cosZ9) 
2fcO80  COS0 


a^  -  cos6(n0^  -  n  cos  9) 

(1) 


sinZ0 


a2  “  n 


cos  9  -  n 
sin^0 


(13) 


(14) 
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Therefore: 


°1 

°2 


cos  6  (n^/  _  -  n  cos  6)  2f 
t  sin^O 

(n„^  cosQ  -  n)  2f 
t  sin^0 


(15) 


Model  epoxy  resin  unidirectional  composites  containing  0.5mm  diameter 
glass  filaments  have  been  cast  inside  glass  tubes  whose  dimensions  are  5mm 
OD  and  3mm  ID.  Prior  to  casting,  all  the  glass  components  were  thoroughly 
cleaned  using  the  following  procedure: 

1 )  Ion-bombardment 

2)  De-oxidation  by  exposing  for  10-15  minutes  to  a  mixture  consisting 

of  170ml  HN03 

30ml  HF 

800ml  distilled  water 

3)  Tap-water  spraying,  10-15  minutes 

4)  Deionised  water  spraying,  2  minutes 

5)  Drying  in  hot  air,  minimum  of  30  minutes 

6)  Cooling  to  room  temperature,  maximum  of  120  minutes 

Small  bundles  of  fibres,  usually  about  10  or  12,  were  gathered  together 
and  inserted  into  the  glass  tube.  In  the  early  stages  of  the 
investigation  no  attention  was  paid  to  the  packing  arrangement,  the  fibres 
usually  adopting  an  irregular  array.  Epoxy  resin  (Araldite  MY753 
manufactured  by  CIBA-GEIGY  (UK)  Ltd)  was  then  driven  into  the  specimen 
until  all  the  air,  including  all  visible  bubbles,  had  been  displaced.  The 
resin  was  gelled  at  room  temperature  for  2  days  and  cured  at  60°C  for  6 
hours  in  accordance  with  the  manufacturer's  recommendations. 

1mm  thick  transverse  sections  were  removed  from  the  specimen  at  intervals 
along  its  length  using  a  high  speed  diamond  impregnated  annular  saw.  In 
order  to  render  the  sections  optically  transparent  and  avoid  corrections 
for  refraction  when  examined  by  oblique  incidence,  it  was  necessary  to 
coat  lightly  the  as-cut  surfaces  with  immersion  oil. 

Figure  17  is  a  photograph  of  one  such  section  examined  with  a  Nikon 
Apophot  180-00-103  polarising  microscope. 

Figure  18  is  an  enlargement  of  a  small  region  taken  from  another  section 
of  the  same  specimen  and  shows  the  complexity  of  the  stress  birefringence 
created  during  cooling  from  the  cure  temperature. 
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The  range  of  stress  birefringence  resolved  in  these  figures  is  large, 
there  being  several  fringe  orders  present.  In  an  attempt  to  reduce  the 
complexity  of  the  stress  pattern  and  to  render  it  amenable  to  comparisons 
with  theoretical  investigations,  specimens  with  parallel  fibres  arranged 
in  a  periodic  array  were  prepared.  Individual  5cm  fibre  lengths  were 
carefully  laid  up  outside  of  the  tube  in  a  rectangular  array  and  their 
ends  locked  into  position  with  quick-setting  epoxy  resin.  The  fibre 
bundle  so  obtained  was  then  cleaned  in  the  normal  way  and  inserted  into 
the  glass  tube  with  the  ends  outside  of  the  tube.  Epoxy  resin  was  then 
driven  into  the  tube  and  allowed  to  gel,  the  fibre  ends  were  cut  off  and 
the  tube  cut  into  two  equal  lengths  which  were  then  placed  in  an  oven  to 
effect  curing.  Using  this  technique  it  was  possible  to  obtain  two  nearly 
identical  specimens,  one  of  which  could  be  retained  as  a  control  specimen. 

Figure  19  shows  photographs  of  1mm  thick  transverse  sections  cut  from 
various  positions  along  a  2cm  length  of  specimen.  The  fibres  were 
arranged  on  a  rectangular  grid  and  photographs  were  taken  using  both 
plane  and  circularly  polarised  light.  Although  stress  differences 
revealed  by  the  isochromatics  vary  non-linearly  along  the  length  of  the 
specimen,  the  general  distribution  of  the  in-plane  stresses  near  the  fibre 
is  consistent  with  finite  element  method  (f.e.m.)  calculations.  It  is 
interesting  to  note  that  large  stress  gradients  surround  some  fibres  and 
fringe  orders  of  at  least  5  can  be  counted. 

Figure  20  shows  a  photograph  of  a  1  mm  thick  transverse  section  cut  from 
the  middle  of  a  specimen  which  exhibited  both  well  defined  fringes  and  a 
large  retardation.  The  image  was  taken  with  plane  polarised  light  incident 
normally  upon  the  specimen. 

Figure  21  shows  a  plot  of  the  hoop  stress  (c^)  vs  distance  from  the  fibre 
for  one  of  the  fibres  seen  in  figure  20.  The  units  of  stress  are  X/Ct. 
With  X  =  546nm,  t  =  1mm,  and  using  C  =  56.6x10  m  /N,  X/Ct  corresponds 
to  a  stress  of  9.65  MN/m  .  The  hoop  principal  stress  at  the  fibre/resin 
interface  is  therefore  of  the  order  of  100  MPa  falling  to  approximately  20 
MPa  one  diameter  from  the  fibre.  The  radial  stress  is  compressive  and  is 
approximately  constant  (6  MPa)  with  distance  from  the  fibre. 

A  specimen  having  similar  geometry  to  that  reported  in  figure  3  was 
immersed  in  distilled  water  at  80°C  for  4  days.  On  removal  from  the 
water,  the  glass  tube  in  which  the  specimen  was  cast  was  observed  to  have 
cracked  at  each  end,  as  shown  schematically  in  figure  22.  Large  stresses 
generated  during  resin  swelling  were  responsible  for  the  cracking.  The 
specimen  was  sectioned  in  the  usual  manner  and  examined  in  the  polarising 
microscope.  The  stress  bi-ref ringence  patterns  are  shown  in  figure  23. 
Very  little  stress  birefringence  was  detected  at  any  section  along  the 
length  between  the  cracked  ends. 
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Whilst  the  general  form  of  the  birefringence  patterns  Indicated  a 
relatively  complex  stress  distribution,  certain  fibres  exhibited  very 
little  asymmetry  in  the  stress  field.  The  birefringence  patterns  observed 
in  figure  20  shows  a  symmetrical  stress  field  with  clearly  defined  radial 
and  hoop  principal  stresses.  The  hoop  stress  (Cqq  *  o^)  is  very  much 
larger  than  the  radial  stress  (c?rr  -  a2).  The  fact  that  arr  does  not  vary 
significantly  with  r  means  that  the  maximum  shear  stress,  (Oqq  “  orr)/2, 
rapidly  increases  near  the  resin/fibre  interface,  i.e.  that  molecular 
orientation,  if  it  occurs,  is  expected  to  be  localised  adjacent  to  the 
fibre.  The  largest  measured  shear  stress,  47  MPa,  was  at  a  distance  of 
some  25  pm  from  the  fibre/resin  interface.  Interpolation  of  the  data  to 
smaller  distances  suggests  that  the  20°C  yield  strength  of  dry  resin 
(~  100  MPa)  would  be  exceeded  at  about  15  pm  from  the  interface. 

The  resin  swelling  that  accommodates  water  uptake  in  composite  materials 
will  depend  on  the  concentration  of  diffused  water  and  is  expected  to  be 
larger  near  to  than  it  is  remote  from  regions  of  water  access.  If 
otherwise  unconstrained,  the  fillet  of  resin  between  three  closely  packed 
fibres  would  endeavour  to  undergo  the  shape  change  sketched  in  figure 
24(b).  Eventually,  the  swelling  in  resin  in  contact  with  the  source  of 
water  molecules  would  saturate  and  the  regions  of  maximum  swelling  will 
move  inwards  as  indicated  in  figure  24(c).  In  the  specimens  studied  here, 
the  resin  was  not  free  to  adopt  these  shape  changes.  The  mechanical 
constraint  offered  by  the  glass  tube  surrounding  the  composite  specimen 
caused  the  resin  swelling  to  generate  a  pressure  of  magnitude  sufficient 
to  burst  the  glass  tube. 


4.2  INTERLAMINAR  SWELLING 

One  source  of  internal  stress  to  be  considered  when  discussing  the  problem 
of  delamination  after  long  in-service  lifetimes  is  inhomogeneous  swelling 
of  interlaminar  resin.  Swelling  of  this  resin  arises  from  the  ingress  of 
water  from  the  outside  and  is  larger  close  to  than  it  is  remote  from  the 
external  surface.  There  is  an  additional  cause  of  inhomogeneity  in  the 
magnitude  of  interlaminar  swelling,  namely  that  due  to  in-plane  anisotropy 
of  the  diffusion  coefficient  for  water  migration;  Kaelble  (ref  -  2nd 
Annual  Army  Composites  Research  Review  1979)  reports  values  of  0.94  to 
10.40  and  0.22  to  0.75  x  10“®cm^sec~^  in  interlaminar  diffusion  parallel 
and  perpendicular  to  the  fibre  direction  in  glass  fibre  reinforced  epoxy. 
(0CF52/Narmco  5205).  In  order  to  investigate  the  interlaminar  stress 
field  attributable  to  water  uptake,  the  following  experiment  has  been 
devised: 

Figure  25  is  a  photograph  of  a  single  laminate  (glass  fibre/epoxy  resin) 
to  which  is  bonded  a  microscope  cover  slip.  By  placing  the  free  surface 
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of  the  cover  slip  in  close  proximity  to  an  optical  flat,  a  pattern  of 
Newton's  rings  is  created  in  the  air  gap  between  the  two.  When  the 
specimen/optical  flat  assembly  is  exposed  to  an  aqueous  environment,  the 
resin  bonding  the  cover  slip  to  the  laminate  swells  and  progressively 
modifies  the  pattern  of  Newton's  rings.  By  superimposing  successive 
photographs  of  the  changing  pattern  on  to  a  photograph  of  the  pattern  seen 
at  the  outset,  it  is  possible  to  create  Moir6  fringes  which  bear  a  1:1 
correlation  to  the  swelling  geometry. 

Figure  26  is  an  example  of  results  so  obtained  for  an  SI /epoxy  laminate 
immersed  in  deionised  water  at  60°C.  The  Moir6  fringes  are  the  fringes 
adjacent  to  the  edges  of  the  square  cover  slip.  Successive  fringes  are 
the  loci  of  displacements  normal  to  the  cover  slip  that  differ  in 
magnitude  by  half  a  wave  length.  The  displacement  fields  parallel  and 
perpendicular  to  the  fibre  determined  after  167  hours'  exposure  are 
summarised  in  figure  27. 

Thin  plate  elasticity  theory  predicts  that  the  4th  differential  of  the 
axial  displacement  (w)  of  a  membrane  that  is  rigidly  supported  at  its  rim 
is  proportional  to  the  pressure  (p)  acting  across  it. 

V4w  -  0  if  p  ■  0 
and  DV^w  «  -p  if  p  t  0 

the  flexural  rigidity 

D*  2  Eh  3 
3(1  -  v2) 

where  2h  *  the  thickness  of  the  cover  slip  and  v  is  Poisson's  ratio. 

Figure  27  shows  a  graphical  representation  of  w  measured  both  parallel  and 
perpendicular  to  the  fibre  direction  after  167  hours'  exposure  to  60°C 
water,  and  the  4th  differential  of  w.  Here  the  maximum  stress  is  3  bars 
tensile.  This  contrasts  with  the  very  much  larger  maximum  stress  of  8.00 
bars  tensile  obtained  with  similar  experiments  on  glass  to  glass  joints 
(section  4.3).  The  maximum  stress  perpendicular  to  the  fibre  direction 
after  the  same  time  of  exposure  to  60°C  water  is  3  bars.  Since  the  origin 
of  the  stress  field  generated  during  water  uptake  is  the  inhomogeneous 
nature  of  the  associated  swelling  (if  the  water  and  hence  the  swelling  was 
uniformly  distributed,  the  resin  would  be  stress-free)  the  fact  that  the 
stresses  between  the  lamina  in  composite  materials  are  apparently  very 
much  smaller  than  those  measured  between  adherends  in  adhesive  joints  is 
attributed  to  experimental  difference  (Redux  312/5  instead  of  American 
Cyanamid  FM1000,  and  640pm  instead  of  150pm  for  the  cover  slip  thickness) 
and  to  a  more  uniform  distribution  of  the  diffused  water.  The  latter  would 
presumably  indicate  a  marked  difference  in  water  diffusion  kinetics. 
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4.3 


INTERFIBRE  SWELLING 


At  the  first  Annual  Army  Composites  Research  Review  (1978),  Professor  D 
Adams  presented  a  finite  element  model  for  the  stress  field  in  a 
unidirectional  composite  that  has  experienced  uniform  differential  thermal 
contraction  between  fibre  and  matrix  materials  and  uniform  swelling  of 
matrix  material  following  water  uptake.  The  assumption  of  uniform 
swelling  is  not  justified;  the  flux  of  diffused  water  received  by  the 
fillet  of  resin  located  between  a  bundle  of  closely  packed  fibres,  is  not 
isotropic  and,  even  for  the  case  of  a  composite  that  is  uniformly 
saturated  with  water,  the  swelling  is  strongly  inhomogeneous  as  indicated 
by  the  sketches  shown  in  figure  24.  Or  special  interest,  because  of  its 
proximity  to  the  resin/fibre  interface,  is  the  Jl  -shaped  film  of  resin 
between  two  adjacent  fibres.  In  the  absence  of  wicking,  the  curved 
surfaces  of  this  film  are  shielded  from  water  by  the  presence  of  the 
fibres,  i.e.  water  uptake  is  restricted  to  diffusion  more  or  less  parallel 
to  these  curved  surfaces.  To  investigate  the  swelling  inhomogeqeity  that 
accompanies  this  pattern  of  water  uptake  and  hence  calculate  the 
inhomogeneous  elastic  field  created  if  (as  is  nearly  true  in  a  real 
composite)  the  resin  is  constrained  not  to  swell,  a  model  specimen 
consisting  of  a  microscope  cover  slip  bonded  with  epoxy  resin  to  a  massive 
and  therefore  rigid  block  of  glass,  has  been  prepared.  As  described  in 
section  4.1  for  the  interlaminar  swelling  specimen,  when  the  free  surface 
of  the  cover  slip  is  brought  into  contact  with  an  optical  flat,  a  pattern 
of  Newton's  rings  is  formed  in  the  gap  trapped  between  them.  When  the 
shape  of  the  gap  changes,  as  happens  when  the  resin  absorbs  water,  swells 
and  displaces  the  cover  slip,  the  pattern  of  Newton's  rings  changes  and 
superposition  of  photographs  taken  before  and  after  the  changes  produces  a 
set  of  Moir'e  fringes  which  faithfully  reveals  the  geometry  of  the  resin 
swelling. 

Figure  28  shows  a  sequence  of  photographs  for  a  square  specimen 
manufactured  using  MY  750  epoxy  resin.  The  cover  slip  was  approximately 
150  /am  thick  and  the  resin  layer  approximately  50 yum  thick.  All  glass 
surfaces  were  thoroughly  cleaned  before  specimen  assembly.  Water  uptake 
was  achieved  by  immersing  the  specimen  in  distilled  water  at  60°C.  From 
the  distribution  of  fringes  it  is  evident  that  the  swelling  is  strongly 
inhomogeneous  and  that  it  faithfully  reflects  the  geometry  of  the  specimen 
for  a  considerable  length  of  time.  It  should  be  noted  that  a  Moir6 
fringe  having  moved  and  been  replaced  by  its  ghbour  represents  an 
increase  in  thickness  due  to  swelling  of  half  a  wavelength,  i.e.  a 
percentage  increase  of  approximately  0.5%.  (The  broad  fringes  running 
across  the  centre  of  the  specimen  and  seen  towards  the  later  stages  of  the 
test  are  a  result  of  the  specimen  tilting  relative  to  the  optical  flat. 
With  improved  temperature  control  and  better  support  for  the  specimen,  the 
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incidence  of  such  fringes  has  been  almost  eliminated).  In  figure  29  the 
position  of  the  first  Moirfe  fringe  is  plotted  as  a  function  of  time. 
Since  this  fringe  corresponds  to  a  relatively  large  water  concentration 
its  migration  is  relatively  slow  and  only  the  early  stages  of  its 
progression  are  recorded.  For  points  near  to  the  rim  and  midway  along  a 
specimen  side,  a  linear  relationship  is  expected  if  the  resin  swelling  is 
governed  by  Fickian  behaviour.  Within  the  limits  of  experimental  error, 
the  relationship  is  linear. 

Using  the  V^w  analysis  outlined  in  section  A. 2,  it  has  been  calculated 
that  the  maximum  stress  developed  in  this  specimen  is  in  excess  of  800 
bars  tensile. 


4.2  FINITE  ELEMENTS 

The  classical  fem  paper  of  Turner,  Clough,  Martin  and  Topp  (1956)  has  been 
developed  and  extended  principally  by  Zienkiewitz  (1977)  and,  for 
application  to  the  microscopic  behaviour  of  fibre  composites,  by  the 
University  of  Wyoming  group.  The  principles  of  the  displacement 
formulation  of  the  fem  for  plane  strain  have  been  applied  to  the  program 
using  triangular  elements  and  linear  shape  functions.  Figure  30  shows  the 
mesh  used  for  developing  the  program.  The  region  bounded  by  nodes  1,  5,  7, 
9,  1  represent  the  fibre  with  the  remaining  mesh  representing  the  epoxy 
resin  matrix.  Figure  31  shows  a  more  realisic  mesh  with  fibre  bounded  by 
1,  19,  35,  1,  where  the  region  19,  36,  52,  35,  19  may  be  treated  as 
interfacial,  again  with  the  remainder  representing  the  matrix.  The  mesh 
of  figure  31  with  236  elements  and  137  nodes  is  generated  logically  in  the 
computer.  This  enables  the  mesh  to  be  modified  for  larger  fibres,  for 
interface  manipulation,  and  for  shape  changes.  This  last  point  is 
regarded  as  particularly  important  since  in  the  form  shown,  i.e.  a  square, 
the  symmetry  may  be  too  high  for  real  composites:  changing  the  right  angle 
of  the  mesh  to  an  arbitrary  value  or  60°  say,  is  easily  achieved.  Both 
mesh  dimensions  are  7  microns  square,  corresponding  to  a  fibre  diameter  of 

9.3  microns  and  a  fibre  volume  ratio  of  39%. 

The  strain  in  an  element  as  a  function  of  the  nodal  displacements  is  given 
by  the  strain  shape  matrix  [B]g,  the  suffix  e  relating  to  a  single 
element.  The  use  of  a  linear  displacement  shape  function  leads  to  a  [B]e 
independent  of  x  and  y,  and  a  constant  elemental  strain.  Linearity  of  the 
dispacement  function  and  the  coupling  of  neighbouring  elements  at  nodal 
points  ensures  continuity  of  displacement  across  element  boundaries,  and 
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compatibility  conditions  are  automatically  fulfilled. 

There  is  a  matrix  [D]g  which  when  multiplied  by  the  strain  vector,  gives 
the  stress  vector  derived  in  terms  of  the  elastic  constants.  The  stress 
system  is  equivalent  to  a  set  of  forces  acting  at  the  nodes  of  the 
element.  The  force  vector  {F>g  may  be  expressed  as 

{F}e  =  [B] g  [D]e  [B]e{a>e 

where  {a>e  is  the  displacement  vector. 

or  {F>e  =  tK]e{a>e 

where  [K]g  is  a  6  x  6  matrix  called  the  elemental  stiffness  matrix. 
Correspomding  to  this  equation  is  a  global  stiffness  equation 

{F>  =  CK]  <a> 

where  the  two  vectors  represent  the  forces  and  displacements  at  every  node 
in  the  mesh.  The  displacement  vector  is  found  by  inversion  of  the 
stiffness  matrix  when  the  loading  condition  is  formulated  in  terms  of 
known  forces  ,  and  unknown  forces  can  be  calculated  from  known 
displacements.  Strains  and  stresses  may  then  be  calculated. 

The  fibre  is  treated  as  purely  elastic  whereas  the  resin  matrix  is 
considered  to  be  stressed  beyond  yield  and  into  the  plastic  region,  i.e. 
resin  plastic  strain  must  be  added  to  resin  elastic  strain.  Thus,  in  the 
program,  loads  are  applied  in  increments  which  follow  the  plastic  stress- 
strain  curve  by  evaluating  a  new  stiffness  matrix  for  each  new  increment. 
Displacements  and  strains  are  calculated  for  each  increment  as  in  the 
elastic  situation.  The  method  is  known  as  the  tangent  modulus  method 
because  the  stress-strain  matrix  is  found  by  using,  instead  of  Young's 
Modulus,  the  tangent  to  the  O' -  C  curve  at  a  given  strain.  The  von  Mises 
yield  criterion  is  assumed  together  with  the  Prandt 1-Reuss  flow  rule. 
Initially  the  whole  load  is  applied  to  solve  for  the  elastic  formulation 
and  if  any  elements  are  stressed  beyond  the  yield  criterion,  all 
quantities  are  reduced  in  proportion  to  bring  them  within  the  elastic 
regime.  The  program  then  proceeds  by  loading  incrementally  where,  at  the 
end  of  each  increment,  the  incremental  strains  are  added  to  the  current 
components  of  strain  and  new  equivalent  stresses  ane  strains  are 
calculated. 

The  stage  has  been  reached  where  the  program  is  functioning  well  for 
elastoplastic  matrix  resin  and  elastic  fibre,  and  the  first  results  for 
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biaxial  tensile  strain  and  pure  shear  are  shown  in  the  figures  32-37. 
Figure  32  shows  the  final  mesh  shape  superimposed  on  the  initial 
undeformed  mesh  for  a  1%  biaxial  tensile  strain.  The  resin  matrix  has  the 
properties  of  an  epoxy  resin  containing  a  water  concentration  of  2%  and  is 
referred  to  in  the  legend  as  low  'plasticity*.  The  crossed  lines  inside 
each  element  indicate  the  orientation  and  relative  magnitudes  of  the 
principal  stresses,  the  program  having  been  written  so  that  the  vector 
lengths  fit  inside  the  smallest  element.  Figure  33  distinguishes  a  third 
phase,  in  that  the  row  of  elements  adjacent  to  the  fibres  have  properties 
of  an  epoxy  resin  containing  a  water  concentration  of  6%,  the  remaining 
matrix  still  containing  2%  water,  and  is  referred  to  as  a  *high 
plasticity'  interface.  Figure  34  is  similar  to  figure  32  except  that  the 
fibre  has  been  enlarged  by  including  an  extra  row  of  elements  in  the 
fibre,  and  demonstrates  the  mesh  flexibility.  Figure  35  is  similar  to 
Figure  33  but  with  the  enlarged  fibre  and  with  the  row  of  elements 
adjacent  to  the  fibre  containing  4%  moisture.  Figures  36  and  37  show  the 
consequences  of  a  1%  extension  in  the  x-direction  combined  with  a  1Z 
compression  in  the  y-direction,  i.e.  the  consequences  of  macroscopic  pure 
shear.  Compressive  principal  stress  vectors  are  barred  at  their  ends. 
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CONCLUSIONS 


1.  By  comparing  optical  retardation  measurements  made  through  parallel 
diameters  in  individual  fibres  in  order  to  detect  interfacial  load 
transfer,  it  has  been  demonstrated  that  the  behaviour  of  epoxy  resin 
composites  during  hot  water  uptake  is  different  for  different  fibre 
chemistries  .  Not  only  are  there  differences  between  'E'  glass  and 
'S'  glass  fibre  composites,  but  also  between  commercial  and 
'contaminant  free'  fibres  of  each  glass  composition. 

2.  S2/SP250  composites,  prepared  with  the  resin  components  mixed  in 
accordance  with  the  manufacturer's  recommendations  for  a 
stoichiometric  mix,  are  not  sufficiently  transparent  for  optical 
retardation  experiments  of  the  kind  reported  in  conclusion  1. 

3.  An  S2/SP250  composite  panel  fabricated  at  AMMRC  is  bounded  by  a  flash 
which  contains  a  low  volume  fraction  of  short  fibres  and  which  since 
it  is  transparent,  is  eminently  suitable  for  optical  retardation 
experiments  to  measuce  times  for  loss  of  load  transfer  during  water 
uptake.  At  80  C  fibres  lying  parallel  to  and  100  m  from  the  external 
surface  lose  their  load  transfer  ability  after  80  hours  exposure  to 
deioniised  water.. 

4.  Using  photoelastic  visualization  of  ultrasound,  it  is  possible  to 
provide  a  ual  display  of  the  effects  of  water  uptake  on  ultrasound 
velocities  in  composite  materials. 

5.  The  swelling  associated  with  water  uptake  by  resin  samples  having 
geometries  representative  of  regions  of  resin  in  real  composites, 
including  composite  laminates,  is  strongly  inhomogeneous.  The 
magnitudes  of  the  corresponding  swelling  stresses  (10's  of  bars)  are 
of  the  same  order  as  those  due  to  differential  contraction  between 
resin  and  glass  during  cooling  from  the  cure  temperature. 

6.  It  is  possible  to  take  account  flow  (due  to  plasticisation  by  water, 
for  example)  when  computer  modelling  the  self-stress  fields  in 
composites. 

7.  Using  measured  differences  in  the  magnitudes  of  radial  and  tangential 
principal  stresses,  it  appears  that  computer  modelling  of  the  true 
(inhomogeneous)  self-stress  fields  of  wet  composites  is  feasible. 
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Figure  1.  Projection  of  the 
E  vector  helix. 


Figure  3.  Effect  of  a  retarder 
on  elliptically 
polarized  light. 
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Figure  4.  Typical  experimental  arrangement. 


Figure  5.  Poincarl  sphere  analysis  of  the  above 


Schematic 


Figure  10  Retardation  in  resin  vs.  radial  distance  from  fibre  centre 
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Figure  16  Photograph  of  a  commercial  'E'  glass  fibre  after 

10200  minutes  immersion  in  water  at  60°C. 


composite  examined  with  the  polarising  microscope. 


Figure  18  Photograph  of  a  small  region  of  a  transverse  section  examined  with 

polariser  and  analyser  crossed. 


Isoclinics  +  isochromatics 


Isochromat  ics 


Figure  19  Photograph  of  transverse  sections  examined  in  both  plane  and  circularlj 
polarised  light.  The  fibres  are  arranged  in  a  periodic  array. 
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Figure  20 

Photograph  of  a  small  region  cf  a  transverse  section  examined  wit 

• 

polariser  and  analyser  crossed. 

Figure  21  Plot  of  the  hoop  tensile  stress  against  radial 

distance  from  the  fibre. 


Fhotograph  of  transverse  sections  examined  in  plane  and  circularly 
polarised  light  .After  immersion  in  water  at  8 0°C  for  4  days. 

The  fibres  are  arranged  in  a  periodic  array. 


tension  compression  ^tension  ;  compression 


Moire  patterns  for  MY  750  resin  immersed  in  water  at  60  C 
showing  the  ingress  of  swelling  during  water  uptake. 
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Figure  32  1%  Biaxial  tensile  strain,  small  fibre, 

low  plasticity  interface. 


Stresses  in  element  1  are:  96.7  M  Pa  and  96.9  M  Pa 


Figure  33  \%  Biaxial  tensile  strain,  small  fibre, 

high  plasticity  interface. 


Stresses  in  element  1  are:  93.  1  M  Pa  and  93.  3  M  P 


Figure  34  1%  Biaxial  tensile  strain,  large  fibre, 

low  plasticity  interface. 


Stresses  in  element  1  are:  145  M  Pa  and  145  M  Pa  . 


Figure  35  1%  Biaxial  tensile  strain,  large  fibre, 

medium  plasticity  interface. 


Stresses  in  element  1  are:  143  M  Pa  and  143  M  Pa. 


Figure  36  1%  Pure  shear  strain,  small  fibre, 

low  plasticity  interface* 

Stresses  in  element  1  are:  47.  7  M  Pa  and  -37.  4  M  P 
Vectors  with  bars  denote  compressive  stresses. 


Figure  37  1%  Pure  shear  strain,  small  fibre, 

high  plasticity  interface. 

Stresses  in  element  1  are  45.  6  M  Pa  and  -36. 
Vectors  with  bars  denote  compressive  stresses. 
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